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Abstract 


The Time-Of-Flight system consisting of plastic scintillation counters plays 
an important role for particle identification in the BESIII experiment at the 
BEPCII double ring ete~ collider. Degradation of the detection efficiency of 
the barrel TOF system has been observed since the start of physical data taking 
and this effect has triggered intensive and systematic studies about aging effects 
of the detector. The aging rates of the attenuation lengths and relative gains 
are obtained based on the data acquired in past several years. This study is es- 
sential for ensuring an extended operation of the barrel TOF system in optimal 
conditions. 

Keywords: 'Time-Of-Flight system, scintillation counter, aging effect, 
attenuation length; 


PACS: 82.80.Rt, 29.40.Mc 


*Supported in part by National Natural Science Foundation of China (U1232201, 
11575225, 11605220, 11675184, 11735014), National Key Basic Research Program of China 
(2015CB856700), the CAS center for Excellence in Particle Physics. 

*Corresponding author 

Email addresses: hhliu@ihep.ac.cn (Huan-Huan Liu), sunss@ihep.ac.cn (Sheng-Sen 
Sun) 


Preprint submitted to Nuclear Instruments and Methods A October 13, 2017 


1. Introduction 


The Beijing Spectrometer (BESIII) detector |1| was designed to study physics 
in the t-charm energy region at the high luminosity Beijing Electron Positron 
Collider (BEPCII) ii . The Time-Of-Flight (TOF) system based on plastic 
scintillation counters plays a key role for particle identification in the BESIII 
experiment, especially for A/a separation; it can also provide fast trigger sig- 
nals. The TOF system consists of a double layer barrel and two single layer end 
caps [1 Al LW In the global BESIII coordinate system, the beam direction is chosen 
as the 2 axis and the polar angle coverage of the barrel TOF is | cos 6| < 0.83, 
while that of the end cap is 0.85 < | cos6| < 0.95. The end cap Time-Of-Flight 
detector was upgraded with multi-gap resistive plate chamber (MRPC) tech- 
nology in summer of 215 ld daH dg In this article, the investigation of the 
aging effects on the performance of the scintillation counters of the barrel TOF 
system of BESIII detector is described. 

The physics data taking of BESIII detector was began in 2009. The detection 
efficiency of barrel TOF system is found to be decreased with time since then, 
and no obvious deterioration of overall time resolution, not “intrinsic” time res- 
olution, is observed. The time dependences of detection efficiency of electrons 
or positrons in Bhabha events for a typical scintillation counter and overall time 
resolution of the barrel TOF system are shown in Fig.[I](a) and (b) respectively. 
The detection efficiency is defined as € = NAt/Ntrack, Ntrack is the number of 
the main drift chamber(MDC) tracks extrapolated to barrel TOF, and naz is 
the number of MDC tracks with their time differences between measured and 
expected times less than a specific value, such as 1 ns, in order to suppress 
the background. The reduction of the number of output signals from the high 
threshold discriminators caused by the aging effect of scintillation counters and 
photomultiplier tubes predominately contributes to the performance degrada- 
tion of the TOF detector. The high voltages (HVs) of the photomultiplier tubes 
(PMTs) of the barrel TOF system were increased twice in 2012 and 2016, the 


reduced output signals were amplified and the efficiency-drops were recovered. 


There are several factors contribute to the overall time resolution of the barrel 
TOF system at BESIII, including the “intrinsic” time resolution, uncertainties 
of time and length of beam bunch, uncertainty of hit position and the expected 
time, time resolution of electronics and uncertainty caused by the time walk [1]. 
Since the contribution of each term was different depends on the performance 
of the detector and accelerator, we could not arbitrarily draw a conclusion that 
the “intrinsic” time resolution did not degrade. The two times increasements of 
HVs enhanced the quantum efficiencies of the PMTs, and were helpful for the 


stability of the number of photoelectrons and the “intrinsic” time resolution. 
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Figure 1: (a) Time dependence of detection efficiency for a typical scintillation counter of the 
barrel TOF system; (b) Time dependence of overall time resolution of the barrel TOF system 
of BESIII. 


Time dependent degradation of detector performance is a phenomenon which 
has been observed for a long time. Serval investigations of aging effects of 
the scintillation counters and PMTs were applied by several particle physics 
experiments, such as the attenuation length aging effect of Belle TOF system |9, 
pa the aging of light yield in the CDF II scintillation counters mE With 
a view to the importance of particle identification and its maintenance over the 
whole lifetime of the experiment, an intensive and systematic study of aging 


effects onto the TOF system is essential for making long-term forecasts. 


2. The barrel Time-Of-Flight system 


The barrel TOF system has two layers of staggered scintillation counters 
mounted on the outer surface of the MDC. Each layer has 88 plastic scintil- 
lation counters hd read out by fine-mesh photomultiplier tubes. The 
barrel counters are Bicron 408 scintillator with a trapezoidal cross-section, and 
the length and thickness are 2300 mm and 50 mm respectively. Two PMTs 
(Hamamatsu R5924) are attached to the two ends of a counter and coupled by 
a 1 mm thick silicone pad (BC634A). The length of the PMT is 50 mm, and 
the total length of PMT assembly is 103 mm, including the base and pream- 
plifier. The TOF readout system consists of preamplifiers mounted in the 
PMT bases, cables, signal time and amplitude measurement circuits and a laser 
calibration system. A dual threshold discriminator scheme, designed to reduce 
noise while maintaining low time walk, is used for time digitization. The out- 
puts from the low threshold discriminators of the two ends are used to start the 
precision time-to-digital conversion (TDC) process and are put into coincidence 
with the high threshold discriminator outputs in order to reduce background 
rates. The signal amplitude measurement circuit is based on the charge-to- 
time conversion (QTC) principle. A more complete description of the Time-Of 
Flight system employed at BESIII, its design and arrangement can be found 
in the articles ne 


the radiation damage, including light transmission, light output, excitation and 


. Since TOF usually operates in a radiation environment, 


emission spectra, of Bicron 408 and other two types of plastic scintillators had 
been experimentally studied in the research and design stage of TOF sys- 
tem at BESIII. The performance of fine-mesh Hamamatsu R5924 PMT had also 
been investigated in a strong magnetic field up to 1T 18 , and an exponential 


dependence of the relative gain on the working voltage was observed. 


3. Attenuation length measurements 


The plastic scintillators exploit the ionization produced by incident charged 


particles to generate optical photons. Organic scintillators, including plastic 


scintillators, do not respond linearly to the ionization density. A semi-empirical 
model by Birks describes the light-output degradation at high ionization den- 
sity : 


dL dE/dx 

de ~ TS Ryd Bae = 
where Z is the luminescence, -% is the luminescence at low specific ionization 
density, dE/dx is the ionisation energy loss in MeV/(g/cm?) and kg is Birks’ 
constant, which is characteristic of the scintillation material. The produced 
photons are then propagated by total internal reflection, refraction and absorp- 
tion in the scintillator until they arrive at the PMTs located at the ends of the 
scintillator and produce a signal pulse. This light attenuation is described as 
an exponential function with an attenuation length A which is defined as the 
length reducing the light signal by a factor of e, and attenuation length is cru- 
cial to achieve a high light collection efficiency for scintillation counters. For the 
barrel TOF, the measured signal pulse height(charge) QTCs at the two ends of 


a scintillator, qı and q2 can be expressed as: 


l/2— z 
q =A: ee “exp (-4 


sin 0 AÀ 
qo l/2+z 
= Áv- . — 2 
Qo = Ay exp ( x i (2) 


where the subscripts 1 and 2 represent the readout channels in the directions 
along the positron beam and the opposite, respectively; Aı and A are the 
relative gains, including the contributions from the light yields of the scintilla- 
tion bars, the quantum efficiencies (QEs), photon-electron collection efficiencies 
(CEs), multiplication gains of PMTs and preamplifier gains for the two readout 
channels respectively; 0 is the polar angle of the incident particle with respect 
to the direction of the positron beam; / is the total length of the scintillator 
bar; z is the extrapolated hit position on the internal surface of the scintillator 
bar along the beam direction based on the particle trajectory 20 ; go is the 
normalized pulse height at z = 0 and 0 = 90°, which is proportional to the 


luminescence of a vertical incident particle. 
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From Eq. [2] it is easily to deduce the following equation 


qı Ay 2z 
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The scattering plot of In(q1/q2) versus hit position z of incident charged particles 
is shown in Fig. P] The attenuation length \ and relative gain ratio A; /A can 
be extracted by means of a linear fitting to In(q1/q2) versus hit position z using 
real data. The nonlinearity behavior close to two ends of the scintillator bar is 


caused by the saturation of QTC measurements and edge effects. 
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Figure 2: Scattering plot of In(qi/g2) versus hit position z. 


The attenuation length dependence as a function of time is expressed as: 
Alt) = Age", (4) 


where Xo is the attenuation length at the beginning of measurements; A(t) is the 
attenuation length at time t; t is the time; a, is the attenuation length aging 
constant for the counter. 

Data samples of electrons and positrons from Bhabha events which are se- 
lected using online event filtering algorithm are employed for the measure- 
ments of attenuation lengths. In the course of several years, measurements of 
attenuation length of each scintillator counter have been performed periodically 
using the method summarized in Eq. B] The plot of time dependence of the 
attenuation length for a typical scintillation counter of the barrel TOF system, 


together with an exponential function from Eq. [4] is shown in Fig. B] The 
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Figure 3: Time dependence of attenuation length for a typical scintillation counter of the 


barrel TOF system. 


attenuation length aging constant is measured individually for each scintilla- 

tion counter. The resulting distribution of attenuation length aging constants 

for 176 scintillation counters of barrel TOF system is shown in Fig. [4] which 

is approximately described by a Gaussian function. The average experimental 
ave 


value of attenuation length aging constants a{’° = 0.045/year, corresponds to 


a degradation annual rate of about 4.4%. 
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Figure 4: Distribution of the individually calculated attenuation length aging constant for 


each scintillation counter together with a Gaussian distribution fitted to the data. 


4. Relative gains 


From 2012 to 2015, the center-of-mass energies of physics data taking were 


2-4.6 GeV, and correspondingly the Bys of electrons and positrons in Bhabha 


events were between 2000 and 4600. Since the relativistic rise of the energy 
loss saturates in this 8y range as described by the “Bethe equations” , the 
differences of mean energy loss per unit length for Bhabha events taken at 
different energies is negligible. That means that the luminescence intensity or 
normalized pulse height qo for Bhabha events at different energies would remain 
constant. In Eq. 2] the relative gains A; and Ao have contributions from the 
light yield of the scintillator bars, efficiencies and multiplication gains of PMTs 
and preamplifier gains of the two readout channels. Presently, it is extremely 
difficult to access the scintillation counters in their positions inside the BESIII 
detector. Therefore it is impossible to extend the investigation of the aging 
effects to a large number of scintillator bars, PMTs and electronics preamplifiers. 
However, since the degradation of the measured signal pulse height is expected 
to be accompanied by the aging effect from the relative gains of each part of 
this sequential chain, the investigation of the overall aging effect is performed 
instead of separate investigations of aging effects of scintillator bars, PMTs and 
electronics channels. 

The products between the relative gains at the two ends of one scintillator 


bar and the normalized pulse height can be derived from Eq. 


Ai-qo = “sin @-exp (4 “)), 


Ee 


(5) 


A2-qo = qg2:sin -exp ( 


Using the same data sample of Bhabha events already used for the study 
of the attenuation length, the distribution between the products of the relative 
gain and the normalized pulse height of a typical electronics channel is shown 
in Fig. [5| together with the best approximation using the Landau function Bi 

The plot of the time dependence of the most probable value (MPV) of the 
product between the relative gain and the normalized pulse height of a typical 
channel is shown in Fig. [6] Since the ideal normalized pulse height qo approx- 
imates to a constant without time dependence, an exponetial behaviour very 


similar to the one already used in Eq. is expected and used to fit the distri- 
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Figure 5: The distribution between the product of relative gain and normalized pulse height 


of a typical electronics channel. 
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Figure 6: The distribution of time dependence of the most probable value of the product 


between the relative gain and normalized pulse height for a typical electronics readout channel. 


70 
60 
50 
40 


Channels 


itr srirtirii Å 
Di -0.05 0 0.05 0.1 0.15 0.2 


Aging constant (/year) 


Figure 7: A Gaussian fit is superimposed to the distribution of the individually calculated 


aging constants for each readout channel. 


bution: 
A(t) = Age" *4* (6) 


where Ao is the relative gain at the beginning of measurements; A(t) is the 
relative gain at time t, t and ay are the time and the relative gain aging constant 
for the electronics readout channel respectively. 

For each electronics readout channel, the relative gain aging constant is ob- 
tained by fitting the data with Eq. [6] The distribution for the aging constants 
of the relative gains for 352 readout channels of barrel TOF system is shown in 
Fig.|7| and this distribution is fitted using a Gaussian function. The average ex- 
perimental value for the aging constant of the relative gains a%”© = 0.032/year, 


corresponds to degradation annual rate of about 3.2%. 


5. Estimation of lifetime of photomultiplier tubes 


The BESIII detector is configured around a 1T superconducting solenoid 
(SSM), considered to be the optimum for precise momentum measurements 
for charged tracks in the 7—charm energy region. Other two superconduction 
quadrupoles (SCQs) are inserted in the conical shaped MDC end caps as close 
as possible to the interaction point. In this region less secondary electrons would 
hit the dynodes and as a consequence reduce the multiplication gain. The fine- 
mesh PMTs R5924 by Hamamastu have a high resistivity against magnetic field 
because of their 19 stages of fine-mesh dynodes, which has many apertures to 
let electrons go through. A fast preamplifier with a gain of 10 is installed in the 
PMT base to boost signals and extend the lifetime of the PMT. 

Since the start of the engineering run of BEPCII in collision mode in 2008 
to 2015, the total running time is about 15869 hours. The average event rate 
is around 300 kHz with the events satisfying the TOF trigger condition, that is 
number of hits in the barrel TOF equal or larger than 1. For each PMT, the 
average hit rate is about 0.04 per event and the most probable value of charge 


distribution is about 12 pC; these numbers are obtained using a random trigger 


10 


data sample. The integrated amount of output charge from individual PMT 
anode is 8.2 C in these 7 years. This is a conservative estimation since small 
signals which could not be over threshold have no record in the raw data. The 
magnitude of integrated charge corresponding to a 7 year time duration is far 


smaller than the limit of lifetime of PMTs, estimated to be around 360 C. 


6. Summary 


Investigation shows that the degradation of detection efficiency of barrel 
TOF system is related to aging effects in the detector, and that no significant 
deterioration in overall time resolution is observed since the start of physics data 
taking of BESIII. Therefore we performed a study on the aging rates of the at- 
tenuation length and relative gain using Bhabha events taken in the period 
of 2012-2015. An approximate estimation of the lifetime of the PMTs shows 
they should be safe under the present irradiation circumstance over the long 
term. These numbers allow us to forecast the detection efficiency and make an 
appropriate plan to ensure the detector operation in optimal conditions. The 
high voltages of the photomultiplier tubes of the barrel TOF system have been 
increased twice since 2009. The relationship between the multiplication gain 
and PMT high voltage could be described using an exponential function and 
the gain will become saturated with enough high HV. Except increasing the 
HV of PMTs, decreasing the thresholds of electronics is another choice for the 
recovery of the detection efficiency. Some experiments about the tuning of the 
threshold are expected to provide a quantitative description of this behaviour. 
The relative gain includes the contributions of light yield of the scintillator bar, 
efficiencies and gain of PMT and preamplifier gain of electronics channel. In this 
study, the aging effect of the relative gain is considered as a cumulative result of 
all these contributions. Individual investigations of aging properties of scintilla- 
tion counters, PMTs and electronics would lead to an improved understanding 
of aging phenomena and enable to extend the operational period, therefore a 


further study is still necessary to obtain reliable long term forecasts, such as 


11 


8-10 years. Thus a detailed proposal for studying the light yield of the scin- 


tillation counters, the quantum efficiency, photon-electron collection efficiency 


and multiplication gain of PMT and preamplifier gain for the readout channel 


is now in preparation. 
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